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ABSTRACT 



We present an affine sl(n+l) algebraic construction of the basic constrained KP hierarchy. 

This hierarchy is analyzed using two approaches, namely linear matrix eigenvalue problem 
on hermitian symmetric space and constrained KP Lax formulation and we show that these 
approaches are equivalent. 

The model is recognized to be the generalized non-linear Schrodinger (GNLS) hierarchy 
and it is used as a building block for a new class of constrained KP hierarchies. These 
constrained KP hierarchies are connected via similarity-Backlund transformations and in- 
terpolate between GNLS and multi-boson KP-Toda hierarchies. Our construction uncovers 
origin of the Toda lattice structure behind the latter hierarchy. 
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1 Introduction 



This paper examines a class of the constrained KP hierarchies and their mutual relations, 
both from the matrix and pseudo-differential point of view. The constrained KP hierarchy 
structure we uncover can be denoted (for n pairs of bosonic fields) by a symbol CKP^ with 
index k ranging from to n and is associated with the Lax operators 
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all satisfying Sato's KP evolution equations. The basic model is the one corresponding to 
k = with the Lax operator M: 
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and the remaining models with k > 1 can be derived from it via the similarity gauge 
transformations at the Lax operator level. Because of the form of the Lax operator (|1.2| ) 
and its abelian first bracket structure we call this CKP^ -* model the n-generalized two-boson 
KP hierarchy. We identify it, at the matrix hierarchy level, with the sl(n + 1) generalized 
non-linear Schrodinger (GNLS) hierarchy ||. Therefore this model can also be denoted as 
GNLS n . We discuss GNLS n model using the structure of the underlying hermitian symmetric 
space and derive its recursion operator as well as its Hamiltonian structure. We also provide 
an affine Lie algebraic foundation for GNLS n hierarchy by fitting it into the AKS sl(n + 1) 
framework ||, generalizing the work of M for the AKNS model. In view of the similarity 
relations connecting this basic CKP^ model and all the other constrained KP models our 
construction uncovers the sl(n + 1) structure behind all the CKP^ hierarchies. 

The other extreme of the sequence CKP^ is the important case of k = n in which we 
recognize the n-boson KP-Toda hierarchy |5|, ||, which recently appeared in a study of Toda 
hierarchies and matrix models |7|, j8|. For one bose pair with n — 1 there is only one model, 
namely the two-boson KP hierarchy || [K]] or CKPj°\ equivalent to the AKNS model. 

In section 2 we recapitulate the Zakharov-Shabat-AKNS (ZS-AKNS) scheme in the frame- 
work of the hermitian symmetric spaces. We work with the algebra Q containing an element 
E = 2/i n • H/a^, which allows the decomposition Q = Ker(adE') © Im(ad-E'). 

We derive a general formula for the recursion and the Hamiltonian operators for a par- 
ticular case where the linear spectral problem XE + A ^ = defining ZS-AKNS scheme 
involves a matrix A e Ker(ad-E). For sl(n + 1) A can be written in the matrix form as 
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In this case we deal with the GNLS n hierarchy. 

We also analyze the ZS-AKNS formalism in case of sl(n+ 1) algebra connected with the 
linear spectral problem with A : 
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which belongs this time to Ker(ad-E) © Im(ad.E'). We give a prescription how to obtain the 
recursion operator in this case, shown in Appendix B (for n — 2) to be equivalent to the 
multi-boson KP-Toda hierarchy recursion operator. 

In section 3 we derive the equations of motion for the GNLS n model using AKS approach 
0. This provides the Lie algebraic foundation for this model, which will turn out to be a 
basic model of the constrained KP hierarchy. In section 4 we prove equivalence of ZS-AKNS 
matrix and KP Lax formulations for the GNLS n hierarchy. Our proof shows that respective 
recursion operators coincide leading to identical flow equations for both hierarchies. 

In section 5 we construct the general CKP^-* models in terms of the GNLS n hierarchy. 
We base the construction on the similarity-Backlund transformation. The structure of this 
transformation ensures that all encountered models satisfy Sato's KP flow equation. In 
particular we obtain a proof for the multi-boson KP-Toda model being the constrained KP 



hierarchy and recover its discrete Toda structure [11 



In Appendix A we collect few Lie algebraic definitions and properties. 



2 ZS-AKNS Scheme and Recursion Operators 

Consider the linear matrix problem: 

Am = d^ (2.1) 

B m ^ = d tr J> 771 = 2,3,... (2.2) 

for 

A = XE + A° with E= 2fia ' H (2.3) 

where \i a is a fundamental weight and a a are simple roots of Q. The element E is used to 
decompose the Lie algebra Q as follows: 

g = JC © M = Ker(adE) © Im(adE) (2.4) 

where the /C = Ker(ad£") has the form K' x u(l) and is spanned by the Cartan subalgebra 
of Q and step operators associated to roots not containing a a . Moreover M. = Im(ad-E) is 
the orthogonal complement of JC. 

In fact, if 2[i a .a/a 2 a = ±1,0, the above decomposition defines an hermitian symmetric 
space Q/K, 0. This special choice of E will play a crucial role in what follows. 
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The compatibility condition for the linear problem ( |2.1| ) leads to the Zakharov-Shabat 
(Z-S) integrability equations 



d m A-dB m + [A, B m ] = ; d = d x d tm = d v 
Inserting the decomposition of A from ( |2.3| ) into ( |2.5| ) we get 

d m A° - dB m + X[E, B m ] + [A , B m \ = 
We search for solutions of (12. tf) of the form 



B n 
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To determine the coefficients of the expansion in 
identity: 

d m A° - Xdm^A + dX(\) + X[E , X{X)\ + [A , X{X)} = 



we first find from ( |2.7| ) the following 

(2.8) 

where X(\) = B rn — \B m -\. By comparing powers of A in ( |2.8| ) we can split X(\) as 

X(X) = O m (l)+Y m (X) ; K m (A)G/C (2.9) 

where O m (l) is an arbitrary element of algebra Q independent of A. Therefore, all dependence 
on A is contained in Y m (X), an element of /C. 

We now consider two distinct cases defined according to whether A is entirely in M. or 
not. 



2.1 A G Im(adE) 

Since now A e M equation ( |2.8| ) yields for the A dependent element Y m in K the condition 
dY m (X) = and hence this component of X(X) can be chosen to vanish. With this choice, 
equation (|2.jj|) becomes a so-called congruence relation: 



B m = XB m ^i + O m (l) = ^2 O m -i{l)X l 



(2.10) 
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Plugging ( 2 . 1 D| ) into ( p.6| ) we find that the coefficients Oj(l) with 1 < j < m — 1 satisfy the 
recursion relation 



80,(1) 
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There is a compact way to rewrite all these recurrence relations for all m. It is provided by 
a single equation: 
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where B is now an infinite series: 
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Since Y m = 0, the equation 
A\ i = 0, 1 



splits into two simple expressions for coefficients of 



d^A = [E,O m (l)} 

d m A° = dO m (l)-[A°,O m (l)} 



(2.14) 
(2.15) 



Equation (|2.14|) can easily be inverted to yield expression for 0^(1). 
From the relation ( |A.5| ) the solution of ( p.!4[ ) is 

0^(l) = [E,d m ^A°] 

From ( |2.15| ) we find now 

90^(1) = [A , 0^(1)} - 0S(1) = d~'[A\ [E 7 d m ^A°] 



(2.16) 



(2.17) 



displaying the non-local character of 0^(1). Plugging now these two results into (|2.15|) we 
find 

d m A° = [E , dd m ^A°] - [A , d~ l [A° , [E , d^A ]]] = TZd m ^A° (2.18) 
where we have defined a recursion operator [S] 



1Z = (d — ad A o d l adjifi^ ad,E 



(2.19) 



Since M. is spanned by 2n step operators associated to roots containing a a only once, i.e. 
a = a a + . . . we denote E±( aa+ „.) = E±, where a — 1, ... ,n. We thus parametrize this model 
by 



A^ = J2(q a E a + + r a E a _ 
With this parametrization we have 



a=l 



0„ 
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where R l sm _ k and R l * are defined in ( |A.8| ). We now find conserved Hamiltonians and the 
first Hamiltonian bracket structure associated to the Z-S equations (|2.5| ) in the Hamiltonian 
form: 

d m A° = { H m , A } (2.22) 

Lemma: d n Tr (A o2 ) = -2dTr (EO n+1 (l)) 

The proof follows from the following calculation: 



d n Tr (A°A°) = 2 Tr ([E, [E, A°]}d n A°) = 2 Tr (e[A°, [d n A°, E]] 
= -2<9Tr = -2«9Tr (EO n+1 (l)) 



(2.23) 



4 



Consider now the Hamiltonian density TCi = YH=iQi r i — Tr(A ) 2 and impose the basic 
relation H, n = d~ l d n 1~Ci for higher Hamiltonians. It follows now from the above Lemma that 
T-L n = — 2Tr (EO n +i(l)), which coincides with the result of Wilson ]nj (see also fl3]). The 
Z-S equations can be obtained from the Hamiltonians defined above and the first bracket 
structure given by 



P 1 6(x - y) 



{ n , Tj } { n , q m } \ = ( -I 
{Qi,r,} {qi,q m } J \I 



5{x - y) 



(2.24) 



Expressing the recursion operator 1Z in terms of the first two bracket structures through 

U = P 2 P{ 1 (2.25) 
we may obtain the closed expression for the second bracket structure of ZS-AKNS hierarchy 



P> 



r s d l r k R 



(d - q s d 1 r k R 
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qsd^qkR 1 ^ 



(2.26) 



This already suggests that the model is bihamiltonian. 

We now consider Q = sl(n + 1) with roots a = oti + c^+i + . . . + aij for some i,j = 
1, . . . , n , E = 2fi n .HVta^, \x n is the n th fundamental weight and H a , a = 1, . . . , n are the 
generators of the Cartan subalgebra. This decomposition generates the symmetric space 
sl(n + l)/sl(n) x (see appendix A for details). 

In matrix notation we have: 



/ 1 



E 



71 + 1 



V 



-11 



J 



This model is defined by A e A4 and is parametrized with fields q a and r a , a 



(2.27) 

.,71 

(2.28) 



which in the matrix form can be written as ( |1.3[ ). The successive flows ( ^.18[ ) related by the 
recursion operator (|2.21 ) are given by 



0„ 



Qi 



n 



(i,l),(j,m) 



(2.29) 
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d, 



n-l 
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From (|2.26| ) an explicit expression for the second bracket is found in this case to be: 

P2-- 



rid 1 r j + Tjd l ri (d - Efc r k d 1 q k ) 5 im - nd 1 q r) 

(d - Efc qkd^rk) S tj - qid~ x Tj qid~ x q m + q m d~ l qi 



(2.30) 
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2.2 A G Im(adE) + Ker(adE) 

The non-zero component A° K in A leads to non-zero terms in equation ( |2.8| ) in K, given by 

d m A K - Xd m ^A° K + dY m (X) + + L4° , Y m (X)) + [A° M , 0^(1)] = (2.31) 

where O^iA) is the component of O n (l) in the subspace M.. 

From equation (|2.31|) we see that Y m (X) = XY m G /C. The A-dependent components in 
subspaces M. and K read respectively: 



-d m ^A° M + [e , o£(i)] + [a^ , y m ] = 

-d m - l A K -dY m + [A° K ,Y m } = o 



while the A independent part leads to 



d m A° M - dO£(l) + [A° M , o£(i)] + K , o£(i)] = 

<9 m A° - SQJJ(1) + , 0*(1)] + [A° , QJJ(1)] 



From ( |2.33| ) we find 



while (|2.32|) gives 
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0^(1) = ad^c^A + ad s ad A o ( £ (d^ad^o V 1 d^d^A 



(2.37) 



It is difficult to present a closed form solution for 0^(1) in general case. However there 
is a simple argument for the existence of an unique solution. Notice that Q2.34j ) provides 
+ 1 algebraic equations (which are not equations of motion) whilst fl2.35| ) provide 



n 



n 



n — 1 such equations making a total of n 2 algebraic equations to solve n 2 = dim JC unknowns 
in O^j. Explicit examples can be worked out case by case (see Appendix B for n — 2). 

Let us consider a particular class of models which are connected to the Toda Lattice 
hierarchy. The models are defined by a Drinfeld-Sokolov linear system of the type 



/ d-X/(n + l) 







d - X/(n + 1) q 2 








\ 





n 



r-2 



d-X/{n+l) q n 

r n d + nX/(n + 1) J 
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with 



A° = J2 qaE {aa) + r a E_ {aa . 



■+a n ) 



= 

(2.38) 
(2.39) 



0=1 



which in the matrix form can be written as in (|1.4j) 
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We now argue that a more general case, with non zero upper triangular elements may 
be reduced to Q1.4Q after suitable redefinition of fields and basis vectors. For this purpose 
consider the generalized DS linear system 



/ d 



qi «i,3 



-1 \ 



d q 2 ■■■ a 2 , n -i 

Qn—l U"n,n— 1 
•■■ d q n 

\r 1 r 2 ■■■ r n d + A J 



(2.40) 



From where we find 
and 



dipn + q n ^n+i = 



(2.41) 
(2.42) 



<9^ n _i + q n -iipn + a nin _iVwi = 
The elimination of if) n +i i n favor of ip n in ( 2.41 ) after reintroducing it in (|2.42|) we get 

^ n _! + qn-^n = (2.43) 



where we have redefined 



and 



?/>n-l = (4>n- 



Q"n,n—1 



q n -i = q n -i + d( 



Qn 

Q"n,n—1 ■ 

q n 



(2.44) 



(2.45) 



We have therefore removed away a n n _i from the linear system. After a succession of such 
redefinitions, we are able to remove all a^j from the linear system showing that ( |2.40[ ) can 
be reduced to (|2.38| ). It can also be shown that the similar construction works in the 
corresponding Lax formalism. 



3 AKS and the Generalized Non Linear Schrodinger 
Equations 

We now derive the GNLS equations from the AKS formalism revealing an afline Lie 
algebraic structure underlying the integrability of the hierarchy studied by Fordy and Kulish 
in connection with the symmetric space sl(n + l)/sl(n) x u(l). This section provides a 
generalization of the work of Flaschka, Newell and Ratiu |3J of the si (2) case. The main 
idea lies in associating an affine Lie algebraic structure (loop algebra Q = Q (g)C[A, A -1 ] of Q) 
to an Hamiltonian structure leading to the same equations of motion as in ( |2.29|) . Here we 
only consider Q = sl(n + 1) leading to GNLS n hierarchy. For the AKS treatment of general 
symmetric spaces see [|14|| . 

The AKS formalism will provide us with a set of Casimir-like Hamiltonians in involution. 
The standard construction is based on the following functions 

fc (X) = iRes(A fc - 1 Tr(X 2 )) ; A = ^ Xj\ j E Q (3.1) 

3 
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where k G Z + . Define now a subalgebra Q of Q as Go = Z),>o^ ® ^ J '- O n the subalgebra 
Qq there is a natural Poisson bracket: 



{f,g}{X) = -(X 



[ir + Vf{X),n+Vg(X)]) ; X e % (3.2) 



where (X\Y) = J2i+j=o ^XjYj and n + is the canonical projection 7r + : £ — ► Go = X)j>o£ <8> 
,V. 

A calculation yields V0 fe (X) = S fc X where S fc is the shift S k : X = £ X^' -> £ X^A^. 
Consequently [V0fc(X) , X] = 0, meaning that 0& is ad-invariant. The main AKS theorem 
ensures now that 0&(X) restricted to Go are in involution: {(pk, <Pk'}(X) = and 
generate commuting Hamiltonian flows. For the parametrization of Go given by Ylj>oXjX ~ J 
with = /ii j) if a + E ( £ a E± a we obtain: 

0fc = E 5AM M + E (3.3) 

i=0 a>0 

where g afe = ^alK~ b , X a fe is the Cartan matrix for sl(n + 1) , E± a are step operators 
associated to the roots ±a. Summation over repeated indices is understood. The Poisson 
brackets fl3.2| ) for h$ and E$ (considered as functions on Go) are 

{h®,h®} = 

{e(a, (3)E^_^ a + (3 is a root 
E:=in a ^ a + /? = (3.4) 
otherwise 

where X Qa = 2a f a , a = I] n a a a and a a , a = 1, . . . , n are the simple roots of sZ(n + 1) PI. 
Hamilton equations in components take the form 

^ = {4> k ,T^} ; T®=h<p,E<p (3.5) 

where the bracket is the Lie- Poisson bracket (|3.4| ). 

The equations of motion are derived from (|3.5|) together with the restrictions for the 
following leading terms 

= o 

= 0, for all roots a 
Eal+...+a a = ?«i a — 1, ... ,n 
E% n+ ^ +aa) = r«, a = l,...,n 

= 0, o,6=l,...,n-l (3.6) 



3 For the sl(n+l) algebra the root system is given in terms of the simple roots as a = ot a + a a+ i~ 
with a — 1, . . . , n and I = 0, 1, . . . , n — 1 (see appendix A for explanation and notation). 
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which subsequently will determine the model. 

From eqn ( |3.5| ) we find the evolution of the step operator E*' to be 



dE U) , /. » 



dt k 2 



-i K n 

l^Ena^E^ + h^E^) 



i=0 a=l 

k 



+ E E {^,(3)E^Et i] + e(-a,m^E^t i] 

i=0 a^/3,a>0 



k n 



+ EE^^^^^ + ^-^^^D (3-7) 



i=0 a=l 



for /3 = X}n a: and 6 1 the Heaviside function. It is easy to see that dE^/d tk = for all k 
once we impose the condition ( |3.6| ) on the right hand side of (577). This shows that = 



v /3 

from fl3.6|) is preserved by the flows. If we now consider k = j = 1 we find after using the 
conditions ( |3.6| ) 

'^ = -±n«hfEf> (3.8) 

m l 0=1 

Consider now a general root not containing a n of the form /3 = J22=l n a a a = a/+ct;+i + ' ■ -+o: s 
with s < n — 1. Since, in this case n n = 0, from the condition (|3.6|), i.e. E^ = we may 



13 

choose hf 1 = for I — 1, . . . , n — 1 allowing therefore a nontrivial "evolution" for E^\ For 
/3 = a n + . . . + a a , and renaming d tl = d x we find 

ES + ... + a a = d x q a (3.9) 
after choosing a normalization = — 1. In a similar manner we have 

E % n+ ... + a a ) = -9,r a (3.10) 
Let us now consider the evolution of h$ given from ( ^.5[ ) by 

^ = EE K ca (EfE^ - E^ES^) (3.11) 

°^A; j=o a >0 

Splitting the sum over the positive roots in terms of its simple root content, i.e. J2 a >o = 
ELlE^ we find 

Qfo(j) k n-ln-s 



X 



V K+-+"Hs) -(o 6 +...+Q 6+s ) -{a b +...+a b+s ) {a b +...+a bA 

UL k j=o s =0 6=1 

x (K b , a + K b+1>a + . . . + K b+S>a ) (3.12) 

It is easy to verify that dh^' /dtk = for a = 1, . . . ,n after imposing the condition (|3.6|) on 
the right hand side of Q3.12| ). Hence the flows preserve the condition = from ( |3.6| ). For 



h\ we find 



E(ftJ5-^+.-+o.) - nE$ b+ ^ +an) ){8 l>b - 5^) = d x (q l+1 r l+ i) - d x ( m ) (3.13) 



9*i 



6=1 
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for I = 1, . . . , n — 1. Due to explicit form of the Cartan matrix K a b, requires an inde- 
pendent calculation yielding 



dtt 



n-l 



,(2) 



(2) 



1=1 



^2d x (q a r a ) - d x {q n r n ) 



(3.14) 



a=l 



By plugging ( |3.9| ), ( |3. 10| ) in (|3.14|) and choosing the integration constants to vanish we get 



h 



(2) _ I qi+m+i - qiri 



I — 1, . . . , n — 1 



(3.15) 



We now determine the evolution of the step operators Ep . From ( |3.7|) we find 



BE, 



(2) 

ai +ai+i + ---+as 

dtx 



(2) 

e(a s+1 + . . . + a n , a t + . . . + a s ) r s+1 E^/ + 



+ e(-(ai + ... + a n ),ai + ... + a s ) qi ^i 2 ( Qs+1+ ... +an) 

Using eqns ( |3.9|) and ( |3. 1C| ) and the relation (|A.2|) implying 

e(a s+ i + . . . + a n , a t + . . . + a s ) = -e(-(ai + . . . + a n ), a t + . . . + a s ) (3.17) 

we may integrate eqn (|3.16|) again choosing the integration constants to vanish. Thus we 
find 

(3.18) 
(3.19) 

(3.20) 



(3.16) 



(2) 

E ai+ai+1+ _ +aa = e(a s+1 + ... + a n ,ai + ... + a s )qi r s+1 
In a similar manner we get 

E m (Qi+ai+i+ ^ +as) = e{a s+1 + ... + a n ,ai + ... + a s )r l q s+l 

Using equation (|3.7|) , (3-9D and (|3.10|) it is not difficult to see that 



0771(2) n 

Pfir, — UI1 '0'n+a n -l + -+aa _ 77.(3) I Q _ „ 



and 



-dlr a 



dE 



dti 



(2) 

-(a„+a„_i+...+a„) 



= 1 



dtx 

Finally, from (|3.7| ) we obtain 

dq a 9E a ^ +a + ___ +aji 



-E 



dU 



(3) 

-(a„+a„_i+...+a a ) 



J -^a a +ot a+1 +...+a n 



2r a £g 6 r 6 (3.21) 



b=i 



and 



dr a _ ^-( ) a a +q a+1 +...+ Qn ) = _ E (3) 



(3.22) 



(3.23) 
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where we made use of relations ( |A.2| ) and between the e's. 

Equations ( |3.21|) and ( |3.20|) yield the following equations of motion after inserting ( |3.22|) 
and ( ggg ) 

d<la r5 2 „ 2„ V* r 
— - C^ga - Zq a qb r h 

dt 2 6=1 

dr a n 
— i = -<9 2 r a + 2r a ^ g b r fc 

for a = 1, . . . , n. These are the GNLS equations ||. 



(3.24) 



4 The n- Generalized Two-boson KP Hierarchy and its 
Equivalence to sl(n + 1) GNLS Hierarchy 

In this section we will establish a connection between GNLS matrix hierarchy for the hermi- 
tian symmetric space sl(n+ 1) and the constrained KP hierarchy. 

The connection is first established between linear systems defining both hierarchies. Re- 
call first the linear problem (2.1), which for A e Ker(ad-E) is parametrized according to 



f d-X/(n+l) 

d-X/(n + l) 







Qi 
Q2 



^2 







r-2 



d-X/(n + l) q n 

r n d + nX/(n + 1) J 



Perform now the phase transformation: 



ip k = exp 



1 



n + 1 



Xdx ipk 



k 



,n + 1 



(4.1) 



(4.2) 



We now see that thanks to the special form of E in A = XE + A , (|4.1|) takes a simple and 
equivalent form: 

( d • • • qi \ ( $i \ 
d ••• g 2 



'02 







(4.3) 



d q n -0, 

\n r 2 ■ ■ • r„ 9 + A y y ^ n+ i J 

The linear problem ( [4.3| ) after elimination of ^ , k = 1, . . . , n takes a form of the scalar 
eigenvalue problem in terms of a single eigenfunction ip n+ i: 



- [d - r kd 1 qk\i>n+i = AV'n+l 



(4.4) 



fc=i 
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This defines the pseudo-differential operator 



which can also be rewritten as 



L n = d - ^2r k d 1 q k 

k=l 



L n = d + ^2ai(d- Si) 1 
i=i 



by a simple substitution: 



n = ate 



-e.-J s > 



i — 1, . . . , n 



(4.5) 



(4.6) 



(4.7) 



The pseudo-differential operator of the form ( f4.5| ) defines a constrained KP hierarchy [1], [15], 

m. 



4.1 On the Constrained KP Hierarchies 

Recall that KP flows with respect to infinite many times t m are defined by: 

d tm L = [(L m ) + , L] {Al 

for the general KP Lax operator: 



L = d + J2ui{{t m })d- 



-l-i 



(4.9) 



i=0 



The subscripts ± mean here that we only take non-negative/negative powers of the differ- 
ential operators d. 



The KP hierarchy (O) allows a straightforward reduction to the so-called /c-th order 
KdV hierarchy by imposing the condition 



L 







(4.10) 



This condition is preserved by the flows in ( [4.8| ) due to 

d tm (L k )_ = [(L m ) + , L k ]_ = [(L™) + , (^)J- (4.11) 
As noted in IT5], [16j condition ( [4.1 0| ) can be made less restrictive by allowing (L fe ) to be 

(L k )_ = -rd- 1 q (4.12) 

since also this condition will be preserved by the KP flows, with the flows of r and q given 
by: 



d tm r = (L m ) + r 



(4.13) 
(4.14) 
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where the upperscript denotes adjoint operation w.r.t. d. In this way one can obtain as a 
special case of ( |4.12|) the two-boson Lax operator: 



L x = d - rdr x q = d + a(d-Sy 1 



(4.15) 



defining a consistent reduction of KP hierarchy with flows defined by d tm L\ = [(L™) + , L±] 
and local bi-Hamiltonian structure || |i~0fl . For the linkage of this restriction of the KP 
hierarchy to the additional symmetries of the KP model see |T7 . 



Obviously in expansion of L\ = d + J2i Uid~ 1 ~ t , the coefficient uq is equal to — rq. This 
can be understood as a symmetry constraint ]1|, [15], IB , which naturally generalizes to 
Uq = — J2k=i r kQk in the case of many bosons. This leads us to the Lax operator in (|4.5| ) 
with flows 



9tmL n 
dt m ri 

dt m qi 



[W) + , 



d~Y, r kd l qu 



k=\ 



k=l / 



We will call the hierarchy defined by the Lax operator L n from 
the n-generalized two-boson KP hierarchy. 



(4.16) 
(4.17) 

(4.18) 



D and flows (^T6|)-(^T8|) 



4.2 Equivalence Between n- Generalized Two-Boson KP Hierarchy 
and sl(n+ 1) GNLS Hierarchy 

We prove now the equivalence between sl(n+l) GNLS hierarchy defined on basis of the linear 
problem ( |4.1|) and the n-generalized two-boson KP hierarchy introduced in the previous 
subsection. Our result is contained in the following 

Proposition. Flows ( \4-lQ )- fy-18i) of n-generalized two-boson KP hierarchy coincide with 
the flows produced by the recursion operator 1Z ^2.2^) of the sl(n + 1) GNLS hierarchy 

We will generalize the proof given in |1| for the n = 1. We proceed by induction. Let 
us first introduce some notation. We parametrize the m-th power of L n (|4.6| ) as 



j<m 



(4.19) 



By defining 



B m = (L^ 1 ) 



(4.20) 



3=0 



the flow equations ( |4.17| )- (|4.18| ) can be rewritten as 



9t m n = B, n 



dt m q t = ~B* m q t 



(4.21) 
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The case m — 1 is obvious. Let us therefore make an appropriate induction assumption 
about m, namely that (|4.21| ) and (|2.29|) agree. Consider now: 



B 



m+l 



d-J2 r kd 1 q k j E p j( m )d J 



k=l 



(4.22) 



n m J 



= dB m + P_i(m) + E E E(-!) V fc (QkPjWf-V 

k=lj=l i=l 

where = d % x f . To calculate the constant term P_i(m) we note that P_i(m) = ResL™ 
which is equal to an Hamiltonian density 7i m and therefore 

/ n \ ra 

P-i(m) = H m = d" 1 ^ - E = - E ^ (ftVt - ^P^) (4.23) 



fc=i 



fc=i 



-EEr 1 (gfcPj-HCrfc)^ - (-l)MP,(m)g fc ) W 
fc=i i=o 



n m ] 



fc=lj=l i=l 

where we have used an identity 



(4.24) 



i=i 

valid for j > 1 and arbitrary A and P. With this information we can now apply P m +i on n: 

n 

Bm+iri = d{B m ri) - E rid' 1 (q k B m r k - r k B* m q k ) (4.25) 
fc=i 

n 

- E 1 (VkBrnri - riB* m q k ) 



k=l 



Writing now the recursion relation of the first section fl2.29|) with induction assumption taken 
into account: 



9m+1 I ft J " U{i ^ \ -B* m q p 



We find from ( ^26|) and (|2~29|) 

n n 

d m+ iTi = -{d{B m ri) - E rid" 1 {q k B m r k - r k B* m q k ) - E r fc^ _1 {qkB m ri - r, ; P^g fc )} (4.27) 

fc=i fc=i 

which agrees with ( f4.25p ( up to a total minus sign). Similarly to ( |4.22j ) we find 

n m j 

-dB* m + p_!(m) + E E E(- 1 ) 5+l ^(^) ( " 1)5 °" l) ^M ( 4 - 28 ) 

k=lj=l i=l 



B 



m+l 
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Applying this on qi we get 



n m J 



B* m+ x<ll = -d(B* m qi) + P-i{m) qi + £E EM)^^)^ (PjHqi)^ (4.29) 

k=lj=l i=l 
n n 

= -d(B* m qi) - qid' 1 (qkB m r k - r k B* m q k ) - ^ qkd~ l {r k B* m q { - qiB m r k ) 

k=l k=l 



following again from identity ( |4.24| ). Again we find agreement (up to a total minus sign) 



with relation (|2.29|) defining recursion operator for the sl(n + 1) GNLS hierarchy. 

5 Constrained KP Hierarchies from GNLS Hierarchy 

In this section we derive a class of CKP models from GNLS hierarchy using the similarity 
transformations. Let us go back to the linear problem defined by A from ( |2.39| ) or ( |1.4j ). 



In matrix notation it is given by Q2.38| ). Performing again the phase transformation (|4.3| ) 



and eliminating ipi, . . . ,ip n from the matrix eigenvalue problem we obtain a scalar eigenvalue 
equation: 

n 1 

- [9 + ^(-l)V,n ^ 1 ^] 4>n+l = \4>n+l (5.1) 
i=l k=i 

This linear problem defines another example of constrained KP hierarchy involving the 
pseudo-differential operators given by: 

n 1 

£>n = 9 + ^(-l)V,n9 _1 ^ (5.2) 

i=l k=i 
n 

= d + Y / a l (d~S l )- 1 ---(d-S 1 )- 1 (5.3) 



i=i 

The coefficients in ( p.2| ) and ( |5.3| ) are related through 

r, = Sl ; q { = -e/^- 1 " 5 *) ; S = ; i = 1, . . . ,n (5.4) 

The pseudo-differential operators of the type shown in (|5.3| ) appear naturally in connection 
with the Toda lattice hierarchy 0. We will refer to the corresponding hierarchy as the 
multi-boson KP-Toda hierarchy. In this section we will explain its status as a constrained 
KP hierarchy and establish its connection with the Toda lattice. 

In |18| it was shown that the first bracket of the multi-boson KP-Toda hierarchy is 
a consistent reduction of the first Poisson structure of the full KP hierarchy. Here we 
will construct successive Miura maps taking the n-generalized two-boson KP system into 
a sequence of constrained KP hierarchies ending with the multi-boson KP-Toda hierarchy. 
This will establish the latter as a multi-hamiltonian reduction of the KP hierarchy. Our 
derivation will reveal the presence of discrete Schlesinger-Backlund symmetry of the multi- 
boson KP-Toda hierarchy. We will also identify the n-generalized two-boson KP hierarchy 
as an abelian structure used in the proof given in Jl8 . 
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5.1 Similarity Transformations 

We start with the Lax operator of the n-generalized two-boson KP hierarchy ( |4.6| ) and define 
a new variable Si = Si + dlnai. Next we perform the following similarity transformation: 

-l 



L?> = (d-S®) L n (d-S?] 

= d + aP^-SPy' + t^id-Sr'id-SPy 1 (5.5) 



i=2 

where we have introduced the redefined coefficients: 



a 



(1) - <?(!)' 



= S[ L) + J> (5-6) 



i=i 



a'p = a[ + di (Si — S[ ) z = 2, . . . , n (5.7) 

If dj = for i = 3, . . . , n we would have already obtained in this way the Lax operator 
of the four-boson KP-Toda hierarchy from (|5.3|) . Otherwise we have to continue to apply 
successively the similarity transformations. In such case, the next step is: 

4 2) - (d-S®) LS> (d-S^ = 9 + ar\d-Si 2 Y + a?(d-S{ 1 Y(d-Si 2 Y 



+ £ af (d - Si)' 1 (d - S^y 1 (d - Si 2) y* (5.8) 

i=3 

= 5 2 + 91n4 1} (5.9) 

4 2) = a^+sf' (5.10) 

n 

4 2) = aW + ^ + ^-^a^+E^ (5.11) 

i=3 

af ] = (d + Si- Si 2) ) af ] 3>i>n (5.12) 

This defines a string of successive similarity transformations. The next step will involve 

similarity transformation (d — S3 ) Lffl (d — S3 J with S3 = S3 + <91na 3 2 ' ) and so on. 
After k steps we arrive at: 

if = a + af» (a - sf;f 1 + af (a - s^)' 1 (a - sf >f 1 + . . . (5.13) 



+flf» (a - si 1 ')- 1 ... (a - sf >)- + £ af> (a - s,)" 1 (9 - si")" ■ ■ • (o- s, 



=fe+l 



Sr=S k + dinar 1 (5-14) 
a? } ee at l) + Sf' (5.15) 

a (fc) = a (*-l) _l ^ _l _ Q( k )\ J*- 1 ) 



(d + s££*> " Si k) ) atf I = 2, . . . , k - 1 (5.16) 

n 

- 4^ + + Sf > - Sf ) afc 1 ) + £ ^ (5-17) 

p=fc+i 

= (d + S p - S? } ) a( fc ~ 1} fc + 1 > p > n (5.18) 
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If n is equal to k + 1 we have obtained in equation ( |5.13| ) the Lax operator L^~^ of the form 
given in ( |5.3| ) and therefore member of the n-boson KP-Toda hierarchy. Since the similarity 
transformations do not change Hamiltonians, the new hierarchy of the Lax operators from 
( |5.3| ) will share the infinite set of involutive Hamiltonians with the n-generalized two-boson 
KP hierarchy. The corresponding Poisson bracket structures are obtained by applying Miura 
maps defined by (CT)-(^T8|). 



At this point we would like to remark on the following ambiguity connected with our 
formalism of successive similarity transformations. Note, that L^ 1 " 1 ^ obtained from ( p,13| ) 



by setting k = n — 1 can be further transformed by an extra similarity transformation without 
changing its form. This is achieved by: 

= d+Y.^ ) {d-sti:?y ■■■{d-s^y (5.i9) 



1=1 



where 



= Sn + dlna^ (5.20) 
a («) = a ^ + S^y (5.21) 

a y = at l) + (o+ styy - sp) aiy i= 2 , . . . , n ^ 

In ( |5.20| )-( |5T21| ) we recognize the Toda lattice structure. Hence the ambiguity encountered in 
associating the multi-boson KP-Toda hierarchy Lax operator to the underlying n-generalized 
two-boson KP hierarchy is an origin of the discrete symmetry of the multi-boson KP-Toda 
hierarchy ||11|| . The discrete symmetry is in this context the similarity map 

L (n-1) _> L (n). 

5.2 Eigenfunctions and Flow Equations for the Constrained KP 
Hierarchies 

Let us start with a technical observation that the similarity transformation, which takes 
£,(*-!) to L^l : 

w = {d - siy L$r* {d - spy 1 (5.23) 

can be equivalently written as 
where we have defined: 

^se/^^aJ^V* (5.25) 

This allows us to find a compact expression for a string of successive similarity transforma- 
tions leading from L n to L^: 



4 fc) = n (*2i^2r l ) L n ff (^-^gr 1 ) (5.26) 
i=k~i v 7 1=0 v 
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It is convenient at this point to introduce a definition of the eigenfunction for the Lax 
operator L. 

Definition. A function $ is called eigenfunction for the Lax operator L satisfying Sato's 
flow equation (\4-8j) if its flows are given by expression: 

9 tm $=(L™) + $ (5.27) 

for the infinite many times t m 

In particular as we have seen in Section 3 the Lax L n (|4.6| ) of the n-generalized two-boson 
KP hierarchy possesses n eigenfunctions $^ = r j = ajexp (/ Si). 

We will now establish a main result of this subsection for the Lax defined in ( |5.13| ). 



Proposition. The Lax operators Lffi ( \5.1c\ ) satisfy Sato's hierarchy equations (U-$ ) and 
possess n — k eigenfunctions given by: 

$W = ag)ef 8 * p=k + l,...,n (5.28) 
Especially we find that = ex P 

(jS { k % 1] ). We will base our induction proof on the 



result of |19j that for given two eigenfunctions $i, $2 of the Lax operator L satisfying Sato's 
hierarchy equation ( |4.8j ) the Lax L = L ^id^ 1 ^^ 1 also satisfies the hierarchy 

equation ( |4.8j ) and the function $1 ^T/ 1 ^) is an eigenfunction of L. 

Start with k = 1. Define = r\ (r^ 1 r^j = (d — (dhxri)) r% for i = 2,...,n. From 
the result of |19[ stated in the previous paragraph and equation ( |5.24|) for k = 1 we find 
that is an eigenfunction of hierarchy of L^\ Furthermore substituting r^'s by a, and Si 
according to ( |5.4j) we find the desired result $^ = ap exp (/ Si). 

Let us now assume that $( fe ) = aWef s f for p — k + 1, . . . , n are indeed eigenfunctions 
belonging to the flow hierarchy of Lffl. Define now 

= (d-(dla^l 1 )^ = (d-(dhia^l 1 )-S w )a^ef^ (5.29) 
= [(d-(dlnai%)-S l+1 + S p )al%I s - = a ( p k+1) el s - (5.30) 



where in ( |5.30|) we have used ( |5.18|) . Since 



L$+» = S&a&Jjr 1 4 fe) O^S" (5-31) 

it follows from induction assumption that L^ +1 ^ satisfies ( |4.8|) . Furthermore since expression 
in ( |5.29p is equal to = $^ ( ^i fe M we have found eigenfunctions belonging to 



p k+1 \ k+l p 

L% +1 > and they are given by (|5 . 2 S ) . This concludes the induction proof. 

As a corollary we find that the Lax operators C n = L^ of the multi-boson KP-Toda 
hierarchy satisfy Sato's hierarchy equations ( |4.8| ) 



Successive gauge transformations of the type given in ([5.24 ) resulting in a Lax structures 



with decreasing number of eigenvalues were also considered in EO] 
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5.3 On the Discrete Schlesinger-Backlund Transformation of the 
Generalized Toda-AKNS model 

In this subsection we will study the discrete Backhand transformations obtained above in the 
Lax formulation within the corresponding matrix formulation of GNLS hierarchy. 

To set the scene let us first connect a "two-boson" Toda system with AKNS matrix 
formulation. The "two-boson" Toda spectral system is: 

(9 -ao(n-l)) = (5.32) 
*n+i + a (n)* n + a x (n)^ n _i = A* n (5.33) 

which becomes in matrix formulation: 

i *i; ) - ( trf ^ 

Under the lattice shift n — > n + 1, \I/ ra goes 

^ n+ i = A^ n - a (n)*„ - ai(n)* n _i = (A - a (n)) ^ n - aa(n)* n _i (5.35) 

Recalling that ao(n) = ao(n — 1) + dlnai(n) we can describe the transition n — > n + 1 in 
form of the matrix equation involving only variables entering ( |5.34j ): 

(0) = (-«iW (A - a (n - 1) - 9 In a x (n)) ) ( *V ) = Tn( *' A) ( V) 

(5.36) 

Compatibility of ( |5.34| ) and ( |5.36| ) yields: 

A n+1 = T n A n T~ l + T n dT- 1 (5.37) 

We will now establish the relation to the AKNS hierarchy. Introduce the following new 
variables: 

and 

q n = -e~ I ao{n - l) ; r n = ai {n)eI ao ^ (5.39) 
Now we can rewrite ( |5.32| )-( [5.33| ) in matrix notation as 



(5.38) 



d q n 
r n d - A 



<>(2) 
ii. 



(5.40) 



If we let <f>$ — > e / X ^ 2 <fin' > we arrive at exactly AKNS equation: 

\ r n d-\\ ){ e -rV2 (i) 



(5.41) 
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Using the Toda chain equations 



da (n) = a\{n + 1) — ai(n) 
dai{n) = ai{n) (a (n) — a (n — 1)) 



and ( |5.35| ) or ( |5.36| ) we find that the shift n — > n — 1 of the Toda lattice results in 



together with 



Qn-l 



X - d In q n q n 
-l/q n 



-r n ql + d 2 q n 



e. 



(<9g n 



r n -i = 



while the shift n — > n + 1 of the Toda lattice results in (equivalently to ( |5.36|) ) 



r" 1 \ / 0i 2 ) 

-r n (A-«91nr n ) J ^ 0« 



together with 



q„ 



q n +i 



r n +i 



6 (2) 



-rlq n + d 2 r n - 



(dr n ) 



(5.42) 
(5.43) 



(5.44) 



(5.45) 



(5.46) 



(5.47) 



These are so-called Schlesinger transformations discussed in ^3|. Here we obtained 

them from a lattice shift of the Toda lattice underlying the AKNS construction. The sim- 
ilarity transformation on the level of corresponding Lax operator captures, as shown in a 
previous subsection, a Toda lattice structure within the continuous constrained KP hierar- 
chy and is fully equivalent to the discrete Schlesinger-Backlund transformation obtained here 
because ( |5.20D and (|5.22| ) correspond to the Toda chain equations (|5.42j ) and (|5.43|) . 

For completeness let us comment on the Schlesinger-Backlund transformations in the 
case of the general Toda hierarchy. We start this time with the spectral equation: 



dV n = V n+1 + ao(n)% 



M 



n—k 



(5.46 



k=l 



The spectral equation (|5.48 ) can be rewritten as a matrix equation, which in components is 
given by: 



f d-a (n-M) -1 

d-a (n-M + l) -1 



\ 









a M {n) 



d-a {n-l) -1 

a M -i{n) . . . a^n) d-X J 



( * 



n—M 



^ i 

* 71—1 



\ 



(5.49) 
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Under the transition n — > n + 1 on the lattice we get: 





-1 \ 






-2 












) 










1 







1 








••• 

V -%(») • • • 



*n-l 

-oi(n) (A-ao(n)) /\*n / 



Introduce now AKNS notation: 



(n) = ^(nje/"^"-^ ; 



J (a (n-i)-a (n— ^ = 1 . M 



(1) 



7 



n— i+1 i 



i = 2, 



M + l 



(5.50) 

(5.51) 
(5.52) 



where for consistency we set ao(n — M — 1) = 0. The relation ( |5.50| ) reads in the AKNS 
variables as: 



( ^ \ I 



l 

o 



o 
i 








••• 

V -r M (n) -r M -i(n) ■ ■ ■ 



r l (n) 
-ri(n) (A — d\nri(n)) J 



(5.53) 

This transformation generalizes the Schlesinger-Backlund transformation from the AKNS 
system to a general AKNS system based on arbitrary Toda lattice. 



A Lie Algebra Preliminaries 

In order to be self contained we recall some basic results on the theory of Lie algebras. We 
first establish the commutation relations in the Chevalley basis, 

{H a ,H b } = 
{H a ,E a } = K aa E a 



e(a, (3)E a+ p a + (3 is a root 
{E a ,E f3 } = { J2 n a=1 n a H a a + P = (A.l) 
otherwise 

where K aa = 2a f a = J2 n bKba, for K ab the Cartan matrix. A root a can be expanded 
in terms of simple roots as a = ^n a a a . The integers l a are defined from the expansion 
_a _ E'""^ e(a,/3) constants are related among themselves by the Jacobi identities and 
the antisymmetry of the bracket. In particular they satisfy 

e( a , /?) = -e(P, a) = e(-/3, -a) (A.2) 
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and 

e(a,/3) = e(/3,7)=e(7,a) (A.3) 

for a + (3 + 7 = 0. 

We choose a special element E in the Lie algebra defined as in terms of fundamental 
weights n a , where 2fi a .ab/a b 2 = Sab as 

E = ^ (A.4) 
or 

The element E decomposes the Lie algebra Q = M. + K, = KeriadE) + Im(adE) where JC is 
composed of all generators of Q commuting with E and is spanned by {H a , E a } with a not 
containing the simple root a a while M. is its orthogonal complement. Hermitean symmetric 
spaces are associated to those Lie algebras Q with roots such that they either do not contain 
a a or contain it only once, i.e. 

5^ = ±1,0 (A.5) 

OL 

a 

for all roots of Q. This fact implies that 

[E , [E , E a }} = E a or (A.6) 

From the Jacobi identities it can be shown that (see e.g. Q/JC is a symmetric space,i.e. 

[JC , JC] e JC, [M , JC] e M, [M , M] e JC, (A.7) 

The Lie algebras satisfying ( |A.5| ) are su(n) , so(2n) , sp(n) , E$, andEj and generate the fol- 
lowing Hermitian symmetric spaces ^gg^j , fgf, ao(10 f° u(1) and ^J^y. 
The curvature tensor associated to these symmetric spaces is defined as 

R s af3 _^Es = [E a , [Ep, E_ 7 ] (A. 8) 

with the hermiticity property 

B sl(3) ZS-AKNS Matrix Model Solution for (1.3) 



Recursion equations for Z-S approach for the case of A^ ' as in (1.4) for n = 2. In this case 
( |2.36| ) gives: 



Y m — —d 1 d m -iqiE ai 



I -dr^d m - x q x 

o o o | (b.i; 
V o o o 



while (|2.371) yields 

("0 O a2 E a2 -\- 0— a2 E_ a2 ± O ai + a2 E ai + a2 ± O — ai — a2 E_ ai __ a2 (B.2) 
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with the coefficients: 



O a2 =d m ^q 2 ; 0_ a2 = -<9 m _ir 2 + r x d 1 d m _ 1 q 1 (B.3) 
O ai+a2 =q 2 d' 1 d m ^iq 1 ; 0_ ai _ Q2 = -<9 m _in (B.4) 



plugging this into algebraic parts of (|2.34| )- (|2.35|) we obtain expression for 

0£(1) = O ai E ai + 0_ ai E_ ai + O hl Ht + O h2 H 2 (B.5) 

with 

1 



O ai = — (dq 2 d 1 d m „ 1 qi + gi9 m _ig 2 ) (B.6) 
q2 v J 

0- ai =d~ l d m -iriq 2 (B.7) 

O fol =d~ l (gi9 _1 9 m _i(g 2 ri) - rig 2 9" 1 9 m _ 1 gi) (B.8) 

O h2 = -d- 1 d m „ 1 q 2 r 2 (B.9) 



After inserting back into the dynamical parts of ( |2.34|) - ([2.35|) the above results lead to the 



recursion operator of the same form as the one for the four-boson KP-Toda hierarchy. 
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